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Abstract 

(N 

The evolution of the entangled muon-neutrino system emerging from charged pion decay 
is explored both in vacuum and in matter. The study is based on a Weisskopf-Wigner type 
wave-packet description. Explicit formulae are derived displaying modulation and attenuation 
of the oscillations due to additional time scales characterising the production process. The case 
of neutrinos disentangled due to the detection of the muon is also considered. 
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' 1 Introduction 

. The analysis of the neutrino oscillation phenomena necessarily includes beyond the propagation, 
also the production and detection stages, since the fmiteness of the space-time distance between 
the starting and the final events is an essential part of the phenomenon. This feature leads to the 
^ ■ smearing of the energy and the momentum of the particles and requires a wave packet description. 

■ Further time scales characterising the production and the detection appear, which modulate and 
^ . damp the fundamental oscillations. These aspects are constantly discussed in the literature of the 
! past decades[U|2j[3llll[5]. Terrestrial observability of decoherence effects was estimated very recently 
O \ by Kayser and Kopp[6j. This is the latest development in the ongoing attempts to extract more 
^ ■ information from the oscillation experiments beyond the oscillation length 0IB]. 

Recent investigations focus on the entanglement aspects of the lepton-neutrino pair propagation 
in different "experimental" situations [9]. Most recently Wu et al. [10] analysed the dependence of the 
density matrix on the detection time of the accompanying lepton when it is detected in addition to the 
propagation time of the neutrino. Their analysis of a two-body decay producing an entangled lepton- 
neutrino pair relies on first order perturbation theory which is strictly valid only for times shorter 
than the lifetime of the mother-particle. This approach seems to be more relevant for long-lived 
sources (for example for beta decay of long-lived radiactive isotopes). Decohering effects reflecting 
the finite width of a Gaussian wave packet were analyzed in Ref . [11] . Decoherence resulting from the 
finite lifetime of the mother-particle was systematically investigated in Refs. [T2| [13] . Both studies 
considered disentangled neutrinos. 

Our goal in this paper is to derive simple expressions governing the coherence of the system in a 
unified treatment which includes the effects of all relevant timescales and relies on an approach which 
makes use of wave packets. In our discussion we make use of two-body decays and the description is 
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valid for times much larger than the source's lifetime. An important process of this type is the decay 
of charged pions tt — > \i + u^, to which we shall mostly refer in this Letter. The treatment will be 
extended also to the propagation through matter where density dependent resonant effects increase 
of the coherence length is observed, which reaches even infinity for some specific value of the lepton 
density. Consequences of muon detection will be carefully explored and compared to the oscillations 
displayed by an entangled neutrino. 

The Letter starts by an overview of the results of the entangled wave-packet treatment of the 
2-flavor muon+neutrino propagation following Ref.[3]. Simple approximate formulae which describe 
the modulation and decoherence of the basic oscillations both in vacuum and matter represent the 
main results of our note. In the second part we repeat the analysis for the case, when the muon is 
detected near the decay tube of the pion. We achieve also in this case a unified picture of the effects 
of finite lifetime and the period of the oscillations observed due to the finite propagation time. 

2 Wave-packet evolution without observing the entangled 
muon 

The state vector of an entangled collinearly outflying muon-neutrino pair in the Weisskopf-Wigner 
approximation reads as follows [3]: 

Y5>(x u ,x^t)) = cos6|mi) + ^ 2 (x u ,x^t) sin9|m 2 ))|/i), 

^(x u , x„ t)=M fdp v l dpMp u + Pfi ) e i(p^+P^)-i(En+^)t. 



E n + E lx -E 1 , + iMj:/2E^ v ' 

where E 2 . = m 2 + pi, E 2 = M 2 + (p u + p^) 2 and the state vector of the muon-type neutrino created 
from the decay is |z/ M ) = cos@|m!) + sinO|m 2 ). The label rrii refers to the mass of the i-th mass 
eigenstate, f n (p) is the momentum profile, T is the energy width of the decaying pion, which will be 
later assumed to have zero momentum on the average. M is an appropriate normalisation constant 
which ensures that / dx v J dx^ip^ 2 = 1 is fulfilled. Eq.([T]) arises from the complete Weisskopf-Wigner 
expression when t » (It should be contrasted with the treatment of Ref.[10] which is valid oly 
for early times t << 

The survival probability of can be studied by projecting ([T]) on the initial (muon + muon- 
neutrino) state. The spatial propagation of the //-type neutrino profile is tracked when the muon is 
not detected, e.g. one integrates the absolute square of the projection over x^. 

Prob{x u ,t) = Jdx^ ^cos 4 6|Vi| 2 + sin 4 e|V 2 | 2 + ^ sin 2 (26^(^2)) . (2) 

Even simpler is the question about the survival of without its localisation, e.g. integrating also 
over x v . The oscillating part of this probability is proportional to the real part of the following 
integral: 



I12 = / dx v J dx^2 = J dp, J dp^\f{p u +p^)\ 2 e l{E ^- E ^ )t 



x i ($) 

(E Vl +E„-E W + iM n T/2E 7T ) (E U2 +E^-E n - iM n T/2E T ) " v ; 
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Another way to recognize its interest is to realise that it determines the off-diagonal element of the 
flavor density matrix in the mass eigenbasis, and characterizes this way the coherence of the two-state 
system: 

pS?W = ^i 2 sin(2G). (4) 

Oscillation phenomena actually are governed by this component. The experimentally observed os- 
cillation in the diagonal matrix elements in the flavor or interaction eigenbase (p^) is just the 
consequence of the orthogonal transformation connecting the two bases. 

A simple and numerically accurate approximate representation of J12 can be derived following [3] 
by expanding all factors of the integrand of ([3]) to linear order around some sort of "mean" momenta 
(denoted by capital letters): P U ,P^, and P n . They are determined from the requirement of energy 
and momentum conservation upon neglecting the mass of the neutrino. After the expansion one 
demonstrates with an appropriate change of variables that the two-variable integral factorizes and 
one integration even can be performed via Cauchy's theorem. The final result is of the following 
form: 



.Am 2 A 2tt 



(12 = exp \ i~2p-t J " „. Am2 eZ F W ex P 




2P V Aur 

F(t) = J rfp|/(p)| 2 exp (- ^p 2 2 ^ ^ (p-Pjt). (5) 

Here Am 2 is the squared mass difference of the two neutrinos, and v = P n /E n , = P M /(P 2 + m 2 ) 1 / 2 . 
The second factor in this expression suppresses the amplitude of the oscillations only when the lifetime 
of the mother-particle is much larger than the oscillation time, which is not the case of the pion decay. 
Clearly, the integral F(t) represents the Fourier transform of the wave function profile of the pion, 
and implies the decay of the oscillating term with characteristic size 1/d. The last factor provides 
attenuation of the oscillation with some 300 years characteristic decay time when one makes use of 
realistic parameters [13] . The coherence length emerging from the concurrence of the last two factors 
is determined by max(T,d). Below we shall assume P n = 0, since then the formulae look much 
simpler. 

Modified formulae reflecting medium effects can be written in an analogous form. Starting with a 
Hamilton operator which contains the Fermi interaction of the e-type neutrinos with the constituents 
of the medium [14"] one solves the Schrodinger equation of the density matrix. For the leading effect 
the mass difference in the denominator can be neglected (in both factors we write for the neutrino 
energy its value calculated at P v ). Then for the probability of the flavor flip one finds 

A 2 \ 2 

Am 



-> u e ,t) = (2tcN) 2 J dp v j d Plx [ Am 2™( p J \f(P»+Pn)\ 2 

(6) 



1 - cos[(Am 2 eff (p u )/2p u )t}) sin 2 (29)/2 



|P Py +P M -^-T72| 2 

Here we have introduced the medium and momentum dependent effective squared mass splitting 

Am 2 eff (p) = ^(Am 2 ) 2 + c 2 p 2 - 2cpAm 2 cos(26), c = 2V2NG F (7) 
{Gp is the Fermi constant, N the density of the medium). 
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Figure 1: The starting (t = 0) probability distribution of a rightward propagating neutrino wave- 
packet. The inverse lengths characteristic for the figures from the left to the right are the following 
(measured in MeV): (T,d) = (0.01,0.001), (0.005,0.01), (0.001,0.01), respectively. 



Assuming slow variation of the effective mass splitting over the momentum range of the neutrino, 
the double integral can be evaluated with the same approximate technique as explained for (j3J). It 
leads to the following factorized form: 



/ Am 2 \ 2 
P(vn -> v e ,t) = — — 2 sin 2 (20)(l - Rel matter ), 

( Am 2 eff (P u ) \ ( \Arh 2 \ T \ 

Imatter = 27T exp U t I F matter{t) exp I ~^pT \ _ ~ t J ' ( 8 ) 

The absence of the second factor of the analogous expression (J^J tells that here we discuss the case 
when the lifetime of the pion is much shorter than the oscillation period. F matt er has the same form 
as in ((Sj) with a substitution Am 2 — > Am 2 , where 

A „ Am 2 ( Am 2 - cP u cos(2Q)) 

A ^ = . o / Try \ -■ 9 

Am 2 e// (P,) 

As a function of c (e.g. the matter density) Am 2 is nonmonotonic and with realistic parameters for 
N 2 x 10 2 iV4/cm 3 it reaches zero (Na is the Avogadro number). It is worthwhile to emphasize 
that the quantities influencing the variation of the oscillation time and of the coherence length in 
matter are different. 

Although with the integration over x v information was lost on the location of the neutrino, 
one might estimate it intuitively with the identification x u ~ t. More detailed information can be 
obtained by evaluating (J2J), which we are going to discuss next. 

In order to ease the triple momentum integral in its last term occuring after the x M integration, 
we have assumed that /(p) is just a window function of size 2d. The result still reflects all generic 
features not depending on the specific form of the pion profile. In this specific case by careful 
mapping of the range of variation the integration can be performed analytically. Denoting its 
(not-quite-transparent) result by G(p v , q v ) one reduces the integral to the form 



I(x v , t) = 2irN / / dp v dq v G(p v ,q v )exp[i(q v —p v )x v -i(E vl (q v )—Eri(p v ))t]. (10) 

J J \p v — qv\<2d 

The remaining two integrations were performed numerically. The shape of the profile of the muon- 
type neutrino signal is obtained by adding to this term also the diagonal contributions. The resulting 
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probability distribution depends on the relation of T and d, as can be seen from Fig. [TJ The shape 
propagates with oscillating amplitude and the forms displayed in the Figure will reappear periodically 
with some changes reflecting decoherence. 

When d « Y the profile is symmetric to some central point which can be identified with the 
actual position of the neutrino. When one goes over to the opposite situation where d » T, the 
time-profile becomes asymmetric. Its forward slope becomes very steep of size ~ l/<i, while the size 
of the back-tail is determined by r -1 . Unfortunately the resolution of the neutrino detectors will not 
allow for a long time the detailed mapping of the neutrino signal. Still this part of our analysis shows 
that the detection rate might be influenced by the way the wave function of the detector particle 
overlaps with the neutrino profile. 



3 Effects of muon detection on the neutrino propagation 

In realistic experiments the detection of the accompanying muons usually provides the estimate for 
the neutrino flux generated in the decay tube of the pions. This disentangling event has an impact 
on the observations of the neutrino oscillations which we shall discuss below. 

Our crude picture of muon detection assumes that the detection consists of measuring the momen- 
tum of the muon by a spatially extended detector. Let its measured value be Q and the measurement 
happens at time t^. This disentangling event produces the following one-particle wave function for 
the neutrino: 

Wfi) — 0i cos 0|?7Ji) + 02 sin 0|m 2 ), 

M**,t) = *i{Q)N I dp ,J^ +Q Ji J ^~ E ^e iE ^. (11) 

J E V)i + Efj,{Q) - + iT/2 

Here rii(Q) is an appropriate normalisation factor for the one-particle wave function of the neutrino. 
The exponent of the last factor is constant, therefore can be omitted. 

We wish to compare the oscillating part of the survival probability in this case to (EJ), therefore 
we integrate over x u : 

I^ s (t) = 27m 1 (Q)n 2 (Q)iV 2 

. \f(Pu + Q)\ 2 e^- E ^ 

Vv (E vl + E^Q) -E^- iT/2) (E u2 + E^Q) - E n + iT/2) ' 1 } 

This integral can be estimated analytically in two limiting cases: d « V and d » V. When 
d « T, the momentum uncertainty is smaller therefore p u varies around P v = —Q. The variation of 
the numerator is faster, therefore the denominator can be taken out as a constant from the integration. 
The requirement of the normalisation of If™ to unity when Am 2 = leads to the cancellation of 
these factors with n\ni. In the opposite limiting case d » F the uncertainty of the energy is small 



therefore the "mean" value of p u is calculated from the energy conservation: P v = M w — \Jrri 2 , + Q 2 . 
The remaining integral is performed via Cauchy's theorem. The results of the integration in the two 
limiting cases are as follows: 

I* id « V) « 27rexp (~i^tj J dp\f(p)\ 2 exp [~^Pt\ , 

7- , . / Am 2 \ 1 / Am 2 r \ . . 

#(«« » F, „ exp {-,—t) ^exp {-—t) . (13) 
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It follows that the coherence length in case of the disentangled neutrino is determined by min(T, d). 
One can easily find that physically this corresponds to the case when the coherence is suppressed by 
the vanishing overlap of the different mass-components of the neutrino wave function. Its character- 
istic time is determined by the difference of the group velocities of the two components. 
In conclusion we summarize the findings of our note: 

• In case of the entangled muon+neutrino pair if V and d are of different order of magnitude 
then the coherence time/length is always shorter than the separation time of the two eigen- 
mass modes. In the disentangled case, however, this latter determines the coherence of the 
oscillations. 

• The dependence of the coherence length on the muon momentum and velocity is slightly dif- 
ferent in the entangled and the disentangled cases (cf. ([5]) and ({TBI ). 

In an experimental situation only some fraction of the muons is detected because of the finite ef- 
ficiency of the muon-detector. Therefore the neutrino beam realistically is a mixture of entangled 
and disentangled neutrinos, which will necessitate some refinement of the analysis when the experi- 
ments become sensitive to such details. The measurement of the coherence length of the oscillations 
might extend the conventional oscillation paradigm in the not-too-distant future. Since there is no 
universal characterisation of the decoherence, future measurements of the decoherence length might 
bring informations on the actual wave-packet profile and also on the degree of entanglement of the 
neutrino beam. 
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